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ABSTRACT 

Using data from the IMACS Cluster Building Survey and from nearby galaxy surveys, we examine 
the evolution of the rate of star formation in field galaxies from z — 0.60 to the present. Fitting the 
luminosity function to a standard Schechter form, we find a rapid evolution of Mg consistent with 
that found in other deep surveys; at the present epoch MJj is evolving at the rate of 0.38Gj/r _1 , 
several times faster than the predictions of simple models for the evolution of old, coeval galaxies. 
The evolution of the distribution of specific star formation rates (SSFR) is also too rapid to explain 
by such models. We demonstrate that starbursts cannot, even in principle, explain the evolution of 
the SSFR distribution. However, the rapid evolution of both Mg and the SSFR distribution can be 
explained if some fraction of galaxies have star formation rates characterized by both short rise and 
fall times and by an epoch of peak star formation more recent than the majority of galaxies. Although 
galaxies of every stellar mass up to 1.4 x 1O 11 M0 show a range of epochs of peak star formation, the 
fraction of "younger" galaxies falls from about 40% at a mass of 4 x 10 10 M Q to zero at a mass of 
1.4 x 10 11 M Q . The incidence of younger galaxies appears to be insensitive to the density of the local 
environment; but does depend on group membership: relatively isolated galaxies are much more likely 
to be young than are group members. 

Subject headings: galaxies: evolution, star formation, cosmic evolution 



1. INTRODUCTION 

The last 10 years have seen an explosion in the num- 
ber of studies of the evolution of star formation in the 
universe. Observations of tens of thousands of high- 
rcdshift objects have been obtained, with redshifts de- 
termined by either spectroscopic or photometric means, 
and star formation rates determined from optical spec- 
troscopy, optical, UV, and infrared photometry, as well 
as X-ray and radio emission. Much of the focus has been 
on the evolution of the star formation rate density (Lilly 
et al. 1996, Madau et al. 1996, see Cucciati et al. 2012 for 
a recent summary) . The evolution of the star formation 
rate density is a fundamental cosmological datum, de- 
scribing the star formation properties of universe as a 
whole, but it contains limited information about the be- 
havior of individual galaxies, whose collective histories 
are responsible for it. However, there have also been a 
large number of studies of the distribution of the star 
formation rates of individual galaxies and their depen- 
dence on galaxy properties, environment, and epoch (e.g. 
Perez-Gonzales et al. 2005, Bauer et al. (2005), Juneau 
et al. 2005, Feulner et al. 2005a,b, Noeske et al. 2007a, 
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Bell et al.2007, Scoville et al. 2007, Karim et al.2011). 
Although there is considerable disagreement about the 
shape of the star formation rate (SFR) distribution at 
a given galaxy mass and epoch, the work of the last 
decade has produced broad agreement about some as- 
pects of the star formation history of galaxies since at 
least z = 1.0. At every epoch the specific star forma- 
tion rate (SSFR — SFR/M sta rs) decreases with galactic 
mass, implying that star formation has a more extended 
history in low than in high mass galaxies. At every mass, 
the SFR and SSFR decline with cosmic epoch at a rate 
which is approximately the same for all masses, suggest- 
ing that the two dependancies are roughly separable, i.e. 
SFR{M gal ,t uni ) = f{M ga i) x f(t uni ). 

Much of the emphasis in recent papers has been on 
the history of mass production: how the mass of stars 
grows with time in particular kinds of galaxies and en- 
vironments, leading to galaxies of the forms and masses 
we see today (e.g. Feulner et al. 2005b, Bell et al. 2007, 
Walcher et al. 2008). However, to understand what is 
responsible for the mass production, we need to under- 
stand the processes which drive the evolution of the star 
formation rates of individual galaxies. One does not need 
to go to cosmological depths to discover that, on aver- 
age, massive galaxies formed the bulk of their stars more 
rapidly than less massive ones; the correlation of spectral 
type with star formation history and the (rough) correla- 
tion of spectral type with galactic mass have been known 
for decades. A recent determination using Sloan data 
has been done by Jimenez et al. (2005). This shift with 
time of star formation from more to less massive galax- 
ies is quite often (and quite confusingly) called "down- 
sizing" . However, this term is better applied as Cowie 
et al. (1996) originally intended, to the purported delay 
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in the formation epoch of less massive galaxies or, more 
precisely, to the delay in the epoch at which star forma- 
tion peaked. These two are quite different phenomena, 
which Neistein et al. (2006), distinguish as archaeological 
downsizing, for the mass dependence of star formation 
timescales, and downsizing in time for the mass depen- 
dence of the epoch of maximum star formation. 

Identifying a delay in galaxy formation, whether mass- 
dependent or not, does require observations at high red- 
shift. In addition to Cowie et al. , claims to see such 
an effect have been made by Feulner et al. (2005), and 
by Noeske et al. (2007b). In all 3 studies the claimed 
discovery rests on the existence of galaxies with specfic 
star formation rates too high to be sustained for the age 
of the universe at the epoch of observation. However, 
very high star formation rates, and very high values of 
SSFR, can be produced for short periods by starbursts, 
and there have been multiple claims that it is starbursts 
rather than delayed galaxy formation which are responsi- 
ble for the observed phenomenon. In this paper we shall 
reexamine these issues, using a sample of field galaxies 
from the IMACS Cluster Building Survey (ICBS- Oemler 
et al. 2012, herafter Paper I). The ICBS is a spectroscopic 
survey of regions surrounding rich intermediate-redshift 
clusters performed using the IMACS wide-held spectro- 
graph on the Baade Telescope. Although targeted to 
clusters, the 30' field of IMACS guaranteed that a large 
fraction of the objects observed would be foreground and 
background galaxies, unrelated to the target clusters, 
and constituting an unbiased sample of the general field 
population. We shall show that starbursts cannot be re- 
sponsible for the observed evolution in SSFR, and that 
the epoch of peak star formation of a substantial fraction 
of massive galaxies was delayed by at least several billion 
years after that of the oldest galaxies. 

The organization of this paper is as follows. In sec- 
tion 2 we describe the data sets we have used, which in- 
clude IMACS spectroscopy, direct imaging, and (in two 
of the four fields) Spitzer 24 /jm photometry. We also 
use nearby surveys to construct local galaxy samples to 
provide zero points for the redshift evolution seen in the 
ICBS. In Section 3 we study the evolution of the galaxy 
populations, including luminosity functions and the evo- 
lution of the specific star formation rate. In Section 4 
we discuss the results, but delay until a following pa- 
per (Gladders et al. 2012, Paper IV) a full analysis of 
the history of star formation in galaxies. Throughout 
this and following papers we shall assume cosmological 
parameters of H Q — 7lkms~ 1 Mpc~ 1 , £l m atter = 0.27, 
VLtot = 1.0. 

2. DATA AND MODELS 

2.1. Data Sets 

We use a number of data sets in the following analy- 
sis. This is intended to be a study of the evolution of 
the general field, by which we mean the totality of the 
galaxy content of the universe, unbiased by group or clus- 
ter membership. Our primary distant sample is drawn 
from the ICBS (q.v. Paper I). The ICBS survey fields 
were selected because the Red-Sequence Cluster Survey 
technique (RCS- Gladders & Yee 2000) which we used 
indicated the presence of a very rich cluster at a suit- 
able redshift (0.30 < z < 0.60); therefore, these fields do 



not represent an unbiased sample of the universe. We 
correct this by removing galaxies in a redshift interval 
about each of the one or two clusters in each field whose 
presence we judge contributed to the RCS detection. We 
do not remove the contents of any other groups or clus- 
ters sufficiently poor, or sufficiently different in redshift 
so that they did not contribute to the RCS signal. The 
results should be a reasonable, even if not perfect, ap- 
proximation to a fair general field sample. This selection 
is described in more detail in Paper I. 

The ICBS covers 4 fields, each about 30' in diameter, 
near the celestial equator. In all fields we have deep mul- 
ticolor imaging in cither the SDSS or Johnson system, 
and very deep imaging in the Gunn r band. In two of 
the fields we have Spitzer 24 /iin photometry of virtually 
the entire field to the confusion-limited depth. The en- 
tire spectroscopic sample contains 6008 galaxies, of which 
4578 comprise our field sample. About 80% of these are 
in the redshift range 0.25 < z < 0.75. We have medium 
resolution (~ 5A) spectra of most of the sample, but for 
1140 objects (946 field galaxies) we only have low dis- 
persion prism (LDP) spectra, which yield lower accuracy 
redshifts (a z ~ 0.01) and no optical line strengths. 

Star formation rates are based on, in order of decreas- 
ing preference, (1) optical emission lines plus 24 /jm 
flux, (2) 24 fim flux, (3) Ha flux, (4) H^flux, and (5) 
[O II] flux. New calibrations of all these SFR measures 
are provided in Paper I; although errors in SFR rates 
increase for the less-preferred measures, all are on the 
same systematic scale. Stellar masses are calculated from 
galaxy luminosities and colors using relations of color vs 
mass-to-light ratio derived from stellar population mod- 
els, assuming a Salpeter IMF, by a method similar, but 
not identical to that of Bell & de Jong (2001). More de- 
tails can be found in Paper I. All masses referred to in 
this paper are stellar rather than dynamical masses. 

We construct two ICBS field samples. The 4-Field 
sample is drawn from all 4 ICBS fields, to a depth of 
r = 22.5 and excludes galaxies with only LDP spectra. 
We will use this sample whenever the analysis concerns 
the environment of the galaxies, for which the LDP red- 
shifts are insufficiently precise. Otherwise, we will use 
the IR-ALL sample, which is drawn from the two fields, 
RCS0221 and SDSS0845 that have Spitzer 24 fj,m pho- 
tometry. This sample includes LDP objects, and goes to 
a depth of r — 23.0. In all fields, only about one half 
of all galaxies above the survey limits were observed- the 
selection of objects being driven by slit mask design- and 
a small fraction of those observed did not yield a reliable 
redshift. Therefore, each galaxy has associated with it a 
weighting which reflects the incompleteness of the cata- 
log at the magnitude and location of the galaxy. Details 
of the ICBS data can be found in Paper I. 

We shall also make occasional use of the AEGIS sample 
(Davis et al. 2007). We have converted their SFR's and 
masses to a Salpeter IMF by multiplying by a factor of 2. 
For objects with detections at 24 ^im , the AEGIS group 
have derived SFR's from Spitzer photometry using the 
calibration of Le Floc'h et al. (2005); we have converted 
these values to that appropriate for our own calibration 
of SFR = f(F 2i ) (q.v. Paper I). To convert SFR's based 
on optical lines to our system, we have remeasured a 
sample of the DEEP2 spectra on which AEGIS is based, 
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Table 1 

Galaxy Samples 



Sample z Area (rj°) N gat 



NGP 0.035-0.045 46.1 1898 

PC2MG a 0.03-0.11 39.9 1362 

IR-ALL 0.20-0.60 0.32 1547 

4 Fields 0.20-0.60 0.63 1976 

AEGIS 0.60-1.00 0.50 2097 



in these models is decidedly non-physical, they have been 
the most popular models for star formation in galaxies 
since at least Tinsley (1972). We take the formation 
epoch, tf orm — 0.9Gyr, and take values of t ranging 
from 1 to 10 Gyr, as well as r = oo, i.e. a constant 
star formation rate. The galaxy properties for these evo- 
lutionary histories were calculated using the Burzual & 
Chariot (2003) stellar models. 



a only galaxies of mass M > 4 x 1O 1O M0 

and have derived a conversion between the AEGIS values 
and our own. 

For comparisons with the properties of the present- 
epoch field, we have constructed 2 local samples. The 
NGP sample include SDSS objects with redshifts 0.035 < 
z < 0.045 in the region of the North Galactic Cap within 
ll h < a < 13 h and 15° < 5 < 45°. SFR's arc de- 
termined from Ha , and masses are calculated from the 
observed SDSS colors, using the calibrations in Paper I. 
This volume has an appreciably lower density than the 
cosmic mean, and is deficient in high-mass galaxies. We 
shall use this sample only for analyses concerning the 
group membership of galaxies. For all other analyses 
we use the PM2GC catalog (Calvi, Poggianti, & Vul- 
cani 2011), which covers the volume 0.03 < z < 0.11, 
10 h < a < 14' 1 45"\ -0?28 < 5 < 0?28. (This sample 
is drawn from a very narrow strip, making the deter- 
mination of group membership very uncertain.) Again, 
SFR's were determined using Ha , and masses derived 
from the reported colors. All of the galaxy samples used 
in this paper are summarized in Table 1. The rcdshift 
ranges are those of the subset of the data actually used 
in this analysis, and the number of galaxies is the num- 
ber within the given redshift range. Additionally, the 
PM2GC sample has been limited to galaxies with masses 
M ga i >4x 10 10 M Q . 

Because star formation rates and masses of all samples, 
both local and distant, have been determined using the 
same mix of methods, described in Paper I, there should 
not be any systematic error in derived galaxy proper- 
ties between the various redshift samples. One potential 
source of systematic error remains. Since, as we shall see, 
mean star formation rates increase rapidly with redshift, 
an incorrect slope of the calibrations of star formation 
rates vs optical and/or infrared indicators would result 
in an incorrect absolute rate of evolution of star forma- 
tion rate. All of our star formation measures are tied to 
the calibration of SFR vs24 /im flux plus optical emis- 
sion lines of Calzetti et al. (2010). Since this calibration 
is tied to the very reliable extinction-corrected hydrogen 
recombination line measurements of star formation rates, 
we think it unlikely to be much in error. 

2.2. Galaxy Evolution Models 

We shall compare our results with two sets of simple 
models for the history of star formation in galaxies. Both 
sets assume solar metallicity and a Salpeter initial mass 
function (IMF). The first set have simple exponential 
decay star formation histories: 

SFR ~ e -(*-*/°™)/ T 

Although the instantaneous turn-on of star formation 
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Figure 1. Galaxy evolution models. Left panels- exponential 
decay models. From bottom to top in the upper and lower panels 
(top to bottom in the middle panel), models with timescales r of 
1.0, 1.5, 2.0, 3.0, 5.0, and 7.0 Gyr, and constant star formation. 
Right panels- delayed exponential models, with timescales, from 
top to bottom in upper and middle panels, bottom to top in lower 
panel, r, of 1.0, 2.0, 3.0, 4.0, 5.5, 6.5, and 12 Gyr. Top panels- 
evolution of star formation rate, in arbitrary units; middle panels- 
evolution of stellar mass, in units of the present epoch mass; bottom 
panels- evolution of specific star formation rate, in units of Gyr -1 . 

Our second set of models assumes star formation his- 
tories of the form introduced by Gavazzi et al. (2002): 

SFR ~ ie _ * /2t 

These are usually called these delayed exponential mod- 
els. For these we calculate galaxy models using a dif- 
ferent set of stellar models. These models use the ob- 
served stellar libraries of Jacoby et al. (1984) in the op- 
tical (~ 3400 to 7400 A), and they were extended to 
the ultraviolet and infrared with the theoretical libraries 
of Kurucz (1993, private communication). They include 
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emission lines formed in HII regions, that were calculated 
using the photoionisation code cloudy (Ferland 1996). 
The nebular component was calculated assuming case B 
recombination, electron temperature of 10 4 K and elec- 
tron density of 10 2 cm -3 . The source of ionizing photons 
was assumed to have a radius of 15 pc and a mass of 10 4 
M© . The effect of using two different sets of stellar mod- 
els is negligible compared to the differences due to star 
formation histories. 

The behavior of the two sets of models is presented 
in Figure 1. Although the star formation histories of 
the two sets of models are very different, the mass and 
specific star formation rate histories of the two sets are 
quite similar, sharing common characteristics, which will, 
in fact, be shared by any models with well-behaved (e.g. 
no starbursts) star formation histories. These include 
the facts that, at any epoch, galaxies with high SSFR's 
must have rapidly evolving masses and slowly evolving 
SSFR's, while galaxies with low SSFR's must have slowly 
evolving masses and may have rapidly evolving SSFR's. 
It is also characteristic of these, but not necessarily of 
every well-behaved model, that the SSFR's within each 
set of models diverge with time, and never cross. This 
will be important in our later analysis. 

3. EVOLUTION OF GALAXY POPULATIONS 
3.1. The Luminosity Function 

We determine the Johnson B band luminosity func- 
tion for our data sample, using a maximum-likelihood 
fitting of a Schechter function, as described in Sandage 
et al. (1979). A preliminary fit to the data, in redshift in- 
tervals of 0.1 between z = 0.13 and z = 0.73 yields values 
for a between -0.90 and -1.35, with no obvious trend with 
redshift, suggesting that using a fixed value a = —1.10 
should be suitable for the entire redshift range. Analysis 
of other high-redshift galaxy samples have yielded val- 
ues of a ranging from -1.03 from the zCOSMOS survey 
(Zucca et al. 2009) to -1.30 from the DEEP+COMBO-17 
surveys (Faber et al. 2007), again with no obvious red- 
shift dependance. In Table 2 we present the results of 
fits to our data, with fixed a = —1.10. 

In Figure 2 we compare the results from Table 2 
with other published values at redshifts in the range 
0.02 < z < 0.90, as a function both of redshift and of 
cosmic epoch. Included are those from the CfA Red- 
shift Survey (Marzke, Huchra, & Gcller 1994), Sloan 
Digital Sky Survey (Blanton et al. 2003), 2df survey 
(Norberg etal.2002), DEEP and COMB017 surveys 
(Faber etal.2007), VIMOS survey (Ilbert et al.2005), 
and ZCosmos survey (Zucca et al. 2009). Although there 
is a modest scatter between different determinations pos- 
sibly due to cosmic variance, the overall trend is quite 
clear: Mg has evolved significantly during the last half 
of the age of the universe, and is evolving most rapidly 
at the present epoch. 

The composite data are well fit by the relation 

M* B = -19.53 - 6.077z + 6.8601z 2 - 2.5740z 3 

Which is shown by the dashed line in Figure 2a. The 
dashed line in Figure 2b represents an evolution rate of 
Mg equal to 0.38/Gyr. This is quite rapid, representing 
a decrease in characteristic luminosity of about a factor of 
2 every 2 Gigayears. We can compare this evolution rate 



Table 2 

Evolution of M* B 



z 






0.10-0.20 
0.20-0.30 
0.30-0.40 
0.40-0.50 
0.50-0.60 
0.60-0.70 




-20.44+°,;!° 
-20.8ll°;i°, 
-20.92+°-Q7 
-21 l4+° : ° 8 

Z1 " L -0.08 

-21 20+ 007 
zi.zu_ . 07 

-21 49+°-° 8 


Table 3 

Predicted Evolution Rate of M' B 




Sample 


dM* B /dt 


Observed 




0.38 



exponential 0.14/0.24 

delayed exponential 0.10/0.20 

delayed exponential + 2 Gyr delay 0.26/0.32 

delayed exponential + 5 Gyr delay 0.32/0.44 



with the expectation from the simple models described 
in the previous section. We take the NGP data set, for 
each galaxy determine the appropriate value of t from 
its SSFR, and, from the corresponding model calculate 
dMg/dt. We evolve all galaxies back in time for a fixed 
time interval, fit a Schechter function to the luminosity 
function at the earlier time, and determine dMg/dt. The 
first three rows of Table 3 compares the observed values 
of dMg/dt, for fixed( a = —1.10) and variable values 
of a, with the values predicted by the exponential and 
delayed exponential models. The predicted evolutionary 
rates are significantly smaller than the observed rate. 

We have obtained these results for only two quite sim- 
ple models for galactic evolution, but one can argue that, 
for any well-behaved model (e.g. no bursts, no special 
epoch of observation) there is a limit to the evolution 
rate of Mg, which depends only of the IMF and the 
epoch of galaxy formation. If the optical luminosity of a 
galaxy is dominated by old stars, then, with to ~ 13Gyr 
and a Salpeter IMF, the evolutionary rate is predicted to 
be dMg/dt ~ 0.10/Gyr. One can obtain a higher rate 
only if the luminosity is dominated by young stars, and 
the star formation rate is falling rapidly. However, if the 
SFR is falling fast, then (for well-behaved models) it was 
much higher in the past, so that the total mass of older 
stars must be large compared to the mass of young stars, 
the luminosity will not be dominated by the young stars, 
and the rapid evolution of the SFR will have minimal 
effect on the evolution of the luminosity. One can escape 
this trap only by significantly reducing the age of the 
galaxies, which will reduce the total accumulated num- 
ber and luminosity of the older stars. In the last two 
lines of Table 3 we present the predicted evolutionary 
rates for the delayed exponential models, for two cases 
in which the start of galaxy formation is delayed for 2 
Gyrs and 5 Gyrs. It appears that a decrease of several 
Gyr in the age of the typical galaxy could explain the 
observed rapid luminosity evolution. 
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3.2. Evolution of the Specific Star Formation Rate 

We now turn to the behavior of the specific star for- 
mation rates of galaxies. There is, in the literature, 
some divergence about the shape of the SSFR distri- 
bution at a given mass and epoch. Following Brinch- 
mann et al. (2004), who claim that the SFR vs mass dis- 
tribution of nearby galaxies has a quite narrow distri- 
bution, Bauer et al. (2005), Noeske et al. (2007a), Peng 
et al. (2010), and Rodighiero (2011) claim to see a similar 
"main sequence" in the distributions of high-redshift ob- 
jects. The Brinchmann distribution, as seen in their Fig- 
ure 17, does have a narrow core among low mass galaxies, 
but even there the 2a tails of the distribution extend out 
by ± a factor of 10 2 , and among more massive galaxies 
with Mg a i > 10 10 M Q even the tight core disappears (see 
also Fig. 15 of Salim et al. 2007). The true breadth of the 
SSFR distribution can also be seen in the lower redshift 
bins of Figure 1 of Feulner et al. (2005a). Given the wide 
range in SFR and SSFR at a given value of Mass, and 
given that all samples are necessarily incomplete at some 
level of SFR, the use of a mean value of either quantity 
seems of limited value. Instead, we will use the normal- 
ized cumulative distribution of SSFR, counting from the 
high end. If we include all objects in doing the normal- 
ization, not just those with detected star formation, this 
distribution should be valid down to the value of SSFR 
at which the sample begins to be incomplete. 

Figure 3a presents the cumulative distribution of SSFR 
for galaxies more massive than 4 x 1O 1O M , in 5 redshift 
intervals, covering 0.0 < z < 1.0. The lowest redshift 
group is drawn from the PM2GC sample, and the higher 
redshift groups are drawn from the IR-ALL and AEGIS 
samples. As has been noted many times before, the dis- 
tribution of SSFR moves steadily to higher values with 
increasing redshift. The models summarized in Figure 1 
all predict an increase of SSFR with lookback time, but 
with a rate of change which decreases with SSFR. They 
also predict a rather hard upper limit to the observed 
SSFR at any epoch, unless the SFR is increasing rapidly, 
as in a burst. These are expected to be quite general 



conclusions, for reasons very similar to those given in 
the previous section for the evolution of Mb- The total 
stellar mass of a galaxy is 

M tot = f° SFR{t)fdt 

The quantity / is the fraction of the mass that remains 
in stars, rather than returning to the ISM during stellar 
evolution. The value of / is about 0.7 for a Salpeter 
IMF, and somewhat smaller for some other IMFs. Then, 
at any epoch the observed specific star formation rate 

SSFR{t) = SFR(t)/ J SFR{t)fdt 

= SFR(t)/f(t- t form ) (SFR) 

If the SFR is not increasing with time, then the maxi- 
mum value of SSFR will occur when SFR(t ) = (SFR), 
meaning constant SFR, and will be 

SSFR 

max — » / . \ 
J \t J form) 

For a Salpeter IMF and an early formation epoch 
SSFR mli x will be about 1 x 10 -10 . Furthermore, the 
evolution rate of SSFR max will be 

d(\og(SSFR))/dt = 0.4343— —t 

J [J- £ f orm ) 

For the same parameters, d(\og(S S F R)) / dt will today 
(t = to) be about 0.05/Gyr. Higher evolutionary rates 
are only possible for values of SFR(t ) < (SFR), i.e. for 
galaxies with declining rates of star formation, and with 
SSFR(t ) < (SSFR). 

Figure 3b presents the fraction of galaxies in the sam- 
ples in Figure 3a with values of specific star formation 
higher than the limiting value expected for constant star 
formation rate and a Salpeter IMF. This fraction in- 
creases from about 2% for local galaxy samples to 21% 
by z ~ 0.9 Furthermore, Fig. 3a shows that the shape 
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Figure 3. a- cumulative distribution of specific star formation rate in samples limited by M ga i > 4 X 10 1() . From left to right: PM2GC 
local sample, ICBS 0.2 < z < 0.4, ICBS 0.4 < z < 0.6, AEGIS 0.6 < z < 0.8, AEGIS 0.8 < z < 1.0. b- fraction of galaxies with values of 
SSFR above the limiting value for a constant star formation rate and Salpeter IMF . 



of the distributions of SSFR changes very little, at least 
out to a redshift of 0.6, contrary to the expectation that 
the highest SSFR galaxies should evolve less rapidly than 
those with lower SSFRs. 

3.3. Are Starbursts Responsible for the Rapid Evolution 
of Star Formation? 

Observations of both the luminosity function and the 
specific star formation rate at earlier epochs show that 
these quantities have evolved faster than expectation. 
Several mechanisms have been proposed to explain this 
high evolution rate. The most common suggestion has 
been starbursts. Beginning with Dressier & Gunn (1983) 
and Couch and Sharpies (1987), starbursts have often 
been invoked to drive galaxy evolution, particularly in 
clusters (see also Bekki 1999, Poggianti et al. 1999, Miller 
et al. 2006, Oemler et al. 2009), but also in field galaxies 
(Le Fevre et al. 2000, Bridge et al. 2007, Bell et al. 2005, 
Bauer et al. (2005), Dressier et al. 2009). 

The main dissent has come from the AEGIS team 
(Noeske et al. 2007a), and, at higher redshift, from 
Rodighiero et al. (2011), which have argued that because 
the SFR vs mass relation is narrow and well-defined at 
all redshifts, only a small number of galaxies could be 
undergoing a large burst of star formation at any time. 
Given the breadth of the SFR vs mass distribution, this 
test is not particularly sensitive to smaller bursts. Nev- 
ertheless, one can make a very general argument that, 
although a starburst can substantially change the instan- 
taneous star formation rate in a galaxy, starbursts can- 
not have a significant effect on the distribution of specific 
star formation rates of a population of galaxies. In Fig- 
ure 4 we present the results from several simulations of 
the effect of adding starbursts to a galaxy population. 
We make the extreme assumption that 100% of galaxies 
have starbursts, either with an amplitude of 3 times the 
average SFR and a duty cycle of 1/5, or an amplitude of 
7 times the average SFR and a duty cycle of 1/10. Nei- 
ther produce a significant enhancement of the high end of 
the SSFR distributions. The reason is simple: starbursts 
cannot, by definition, change the time averaged SFR. 
Periods of enhanced SFR must be balanced by periods 
of depressed SFR. For any reasonably shaped SSFR dis- 



tribution, the latter effect offsets almost completely the 
former. 

This is a sufficiently important point that it is worth 
repeating another way. One can contrive any history 
of star formation that one wants, driven by starbursts 
or any other conceivable process. However, the specific 
star formation history is not so easily changed, since it 
is a ratio which depends in its denominator on the entire 
past history of star formation. Unless one is observing 
a special epoch, before which starbursts did not occur, 
adding bursts of star formation to a galaxy's history will 
also add mass and therefore will not by itself change the 
time-averaged specific star formation rate. Stated yet 
another way, the SSFR measures the characteristic (in- 
verse) timescale for star formation, not the star formation 
rate. 

This is not to say, however, that starbursts, if common, 
cannot effect the SSFR values of individual galaxies. The 
red curves in Fig. 4 present, for the two starburst sce- 
narios, the fraction of objects at a given value of SSFR 
that are observed to be in a starburst. Starbursts do not 
have a large net effect of the SSFR distribution because 
the starburst cycle moves objects to lower as well as to 
higher values. However, those with the highest values at 
any time are those currently in a starburst phase. Al- 
though they do not emphasize the fact, Bell et al. (2005) 
find a factor of two higher merger rate among galaxies 
with the highest values of SSFR; and, Bridge (2005) find 
that, at0.2<z<1.3, galaxies with detected 24 /jm flux 
are a factor of 5 more more likely to have close compan- 
ions than those without 24 /im flux. In a forthcoming 
paper (Oemler et al. , in preparation), we shall show that 
a large fraction of ICBS galaxies with very high values 
of SSFR are undergoing starbursts driven by galaxy in- 
teractions. 

3.4. Are Some Galaxies "Young"? 

We have shown that starbursts are an unlikely expla- 
nation for the too-rapid evolution of galaxy luminosi- 
ties and star formation rates. We have argued above 
that lowering the mean age of galaxies can explain 
the fast evolution of Mg, and Cowie et al. (1996) and 
Noeske et al. (2007b) have argued that the same change 




SSFR 



Figure 4. The black lines, with the scale at the left of the graph, 
show the effect of starbursts on the distribution of SSFR in a galaxy 
population. Solid line- Local SSFR; dashed line- population with 
time averaged SFR of the Local population but galaxies have 3:1 
starbursts with a duty cycle of 1/5; dotted line- population with 7:1 
starbursts and 1/10 duty cycle. The red lines, with the right-hand 
scale, show the fraction of objects in the two starburst scenarios, 
which are observed in the starburst phase, as a function of SSFR. 

can explain the rapid evolution of the SSFR (Feulner 
et al. 2005b have raised the same possibility). In this sec- 
tion we will examine this possibility, using the ICBS and 
local data samples. Figure 5a presents the cumulative 
distribution of SSFR for galaxies in the Local data set, 
and the 0.40 < z < 0.60 redshift range in the IR-ALL 
data set, down to a limiting mass M ga i > 4 x 10 10 M Q , to 
which both data sets are essentially complete. We wish 
to compare the evolution shown in Fig 5a to the predic- 
tion of the exponential and delayed exponential models, 
as presented in Fig. 1, both for an early time and for 
reduced ages of formation. 

The fraction of objects with values of SSFR too high 
to explain by early-formation models is one indicator of 
young galaxies, but it is an insensitive test, since even 
quite young galaxies can have low values of SSFR if 
the SFR's fall from their peak values sufficiently rapidly. 
Furthermore, this criterion depends on an absolutely cor- 
rect calibration of star-formation rates, both observed 
and predicted by models. Instead, we will use the evolu- 
tion rate of the SSFR distribution as our main tool. It 
is a characteristic of all of the models in Fig 1 that the 
curves of SSFR(t) diverge with time, and never cross. 
Therefore, given a set of galaxies at one epoch, ordered 
by SSFR, evolution with time will change the values of 
SSFR but will not change the rank order of the galaxies. 
Thus, the galaxies in the Local sample in Fig 5a, at any 
given value of /, should have evolved from the galaxies in 
the IR-ALL sample at the same value of /. Furthermore, 
if either set of models correctly describes the behavior of 
real galaxies, the value of SSFR(t a ) at that value of 
/ will correspond to a unique r model, from which we 
can predict the value of SSFR(t) at the epoch of the 
earlier observations. There are two possible exceptions 
to these statements. Firstly, they ignore the growth of 
mass with time, which will cause some galaxies which are 
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Figure 5. a- The cumulative distribution of specific star forma- 
tion rate in samples limited by M ga i < 4 X 10 10 . left- PM2GC 
sample, right- IR-ALL sample with 0.40 < z < 0.60. b- evolution- 
ary rate of SSFR vs SSFR(t a ). The filled red points correspond 
to the values of f shown as red points on the local population curve 
in the top plot. Open red circles are corrected for the effect of 
mergers, as described in the text. The two grids correspond to 
exponential (solid) and delayed exponential (dashed) galaxy evo- 
lution models with indicated values of tau, and indicated delays of 
the start of star formation (in Gyr) . 

above the mass limit at the present epoch, to be below 
the mass limit, and excluded from the samples, at higher 
rcdshifts. However, we have simulated the effect of this 
on real data samples, and find that it has a negligible 
effect on the SSFR distribution. Secondly, they ignore 
the effect of mergers, which violate the assumption that 
a fixed population of objects evolves from one epoch to 
another. We will return to this issue shortly. 

We define Z(SSFR(t )) = d{log{SSFR))/dt, with 
units of Gyr~ l , and measure its value using the two data 
sets in Fig 5a, at multiple values of / shown by red cir- 
cles. In Figure 5b we compare the measured values of Z 
at these points with the predictions of the exponential 
and delayed exponential models from Fig 1, both for the 
early formation epochs (heavy solid and dashed lines) 
and for various delays in the formation epoch (lighter 
solid and dashed lines). Included are typical la errors, 
calculated assuming random subsets drawn from an over- 



all distribution of SSFR values with a shape like that 
of the observed galaxy samples. The results are incon- 
sistent with a single epoch of formation for all galaxies 
with either set of star formation histories. However, at 
lower values of SSFR, beyond the 25th percentile (count- 
ing from high to low), the results lie close to the delayed 
exponential early-formation models. At higher values 
of SSFR, the data move steadily away from the heavy 
line, suggesting more and more recent epochs of forma- 
tion. Thus, a reasonable conclusion is that most mas- 
sive galaxies (M ga i > 4 x 10 10 M Q ) are coeval and old, 
but that about one quarter, those with values of SSFR 
greater than about 10~ n are "younger", with an epoch 
of formation which moves to later times with increasing 
SSFR. 

To examine the effect of mergers on this result, we 
perform a Monte Carlo simulation, using data from the 
IR-ALL sample in the redshift range 0.20 < z < 0.40. 
We take a merger rate versus redshift and galaxy mass 
as described in Xu et al. (2012) Equation 9. We assume 
that, after the merger, the mass is the sum of the mass 
of the two merging galaxies and that, after the merger is 
complete and any starburst has decayed, the total star 
formation rate is the sum of the pre-merger SFR's of the 
components. Although this is undoubtedly an oversim- 
plification, a more elaborate treatment is not justified by 
our limited knowledge of the phenomenon. With these 
assumptions, the principal effect is a slight steepening of 
the SSFR distribution due to the averaging of the ex- 
tremes of low and high SSFR values as pairs of galaxies 
combine. The effect on Z is shown by the open circles 
in Fig. 5b. Except at very low values of SSFR (in which 
we are least interested) its amplitude is very modest and- 
given our ignorance of the details of the process- we shall 
ignore it in the remaining analysis. 

At this point we need to elaborate on what we mean by 
"delayed formation" . The behavior of real galaxies which 
neither set of models can reproduce is the existence of 
galaxies with both high and rapidly falling SSFR's. Ex- 
ponential models peak at t = tf orm ; those with high 
values of SSFR at later times are those whose SSFR de- 
clines slowly. Delayed exponential models with arbitrar- 
ily high values of r can peak at arbitrarily late times, but, 
like the exponential models, those with high values of 
SSFR at late times are those with slowly varying SSFR. 
It is obvious that any models (other than contrived cases 
seen at special epochs) which peak early must have either 
low values of SSFR or low values of d(SSFR)/dt at late 
times. Only galaxies whose star formation histories peak 
at an age of the universe significantly greater than the 
timescales on which star formation rises and falls before 
and after the peak can reproduce the behavior seen in 
our data. These we call "young" or "delayed formation" 
galaxies. It must be emphasized that, by delayed forma- 
tion, we are only talking about a delayed onset of the 
bulk of star formation. There is no way, with existing 
data, that one can rule out a tail of star formation in all 
galaxies extending to quite early times. Thus, observa- 
tions suggesting that there is some old stellar population 
in most if not all nearby galaxies (Grebel & Gallagher 
2004, Orban et al. 2008 ) does not contradict our picture. 

Cowie et al. (1996) and Noeske et al. (2007) observe 
downsizing: not only are some galaxies young, but their 
age decreases with decreasing mass. To examine this 



effect, we divide our samples by mass. Figure 6 re- 
peats the analysis in Fig. 5, but dividing galaxies into 3 
mass ranges, 4 x 10 10 M o - 7 x 1O 1O M , 7 x 10 10 M o - 
1.4 x 10 n M Q , and M gal > 1.4 x 10 n M o . There is a 
clear divide, at approximately M ga i — 1.4 x lO n M0. 
Galaxies above this limit follow the old delayed ex- 
ponential locus; but at least some galaxies below this 
limit need a range of younger ages. Galaxies in the 
mass range 7 x 1O 1O M — 1.4 x lO n M diverge from 
the old locus below the 30th percentile, galaxies in the 
4 x 10 10 M Q — 7 x 10 10 M Q mass range depart from the 
old locus below the 40th percentile, a very large fraction 
of young galaxies, indeed. Note, however, that at every 
mass less than M ga i = 1.4x 10 11 Mq there is a large range 
of age, with a majority of old, coeval galaxies. Although 
mean age is a function of mass, such an averaging hides 
a more complicated story. 
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Figure 6. a- cumulative distribution of specific star formation 
rate in 3 mass ranges. Left 3 curves- PM2GC local sample; right 3 
curves- IR-ALL sample 0.4 < z < 0.6. Blue curve- 4 X 1O 1O M < 
M < 7 X 10 10 M Q , green curve- 7 X 10 10 M o < M < 1.4 X lO n M , 
red curve- M ga i > 1.4 X 10 1:l Mq b- evolutionary rate of SSFR vs 
SSFR(t ). Colors correspond to mass ranges in Fig 6a; the grid 
of models is the delayed exponential grid from Fig. 5b 
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The expected errors in this case are more significant 
than those shown in Figure 5b. Evaluating the over- 
all significance of the departure of, for example, the red 
and blue curves from these typical errors is not straight- 
forward, because the curves are derived from cumulative 
distributions and the errors of adjacent points are, there- 
fore, at least partially correlated. Instead, we perform 
Monte Carlo simulations on pairs of distributions, drawn 
from the same parent population, and of the size of our 
samples, and ask how often the two resulting Z curves 
depart from each other, at every value of /, by as much 
as the pairs of Z curves in Fig 6b. We find that, for 
example, the chance that the red and blue curves were 
drawn from the same population is only 2 x 10~ 5 , so the 
observed differences are very significant. 
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Figure 7. Same analysis as in Fig. 6, but with galaxies from 
the PM2GC local sample, and the 4-Fields 0.4 < z < 0.6 sample, 
divided into intervals of local density. Blue- lowest third values of 
p/Pmed; Green- middle third values of p/p TOe di Red- highest third 
values of p/p me d 

One might expect that galaxy formation and evolu- 
tion would proceed more rapidly in high-density regions 
of the universe than in low-density regions (e.g. Cen 
2011). Therefore, in Figure 7, we test for dependance 
of galaxy evolution on local density which we calculate 
from the distance to the third nearest neighbor with 
mass Mg a i > 4 x 1O 1O M . We divide the distribution 



of p/ Pmed, where p me d is the median value of p for galax- 
ies in the same redshift interval, into three parts, and 
examine the evolution of SSFR for each third. The 
nearby sample is the PM2GC sample, and the distant 
sample consists of galaxies in the 4Fields sample with 
0.4 < z < 0.6. Because, as we have shown, the evo- 
lution of SSFR depends on galaxy mass, and because 
the distribution of mass is dependent on local density 
(Vulcani et al. 2012a) we weight galaxies in each density 
interval by mass to produce weighted distributions with 
mass which are the same in all intervals. Although the 
behavior of galaxies in the lowest third of the density 
distribution differs a bit from the others, all 3 groups 
seem to display the same wide range in formation epochs, 
and thus provide no evidence that the age distribution 
of galaxy populations depends on local density. 
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Figure 8. Same analysis as in Fig. 6, but with galaxies from the 
NGP local sample, and the 4Fields 0.4 < z < 0.6 sample, divided 
between isolated galaxies (blue), all group members (green), and 
members of groups with 4 or more members (red) 

Figure 8 presents the behavior of the same redshift 
intervals analyzed in Figures 6 and 7, now divided by 
group membership. We have constructed a group cata- 
log, for the redshift interval 0.2 < z < 0.6, from the 4 
Fields sample using the standard Huchra & Geller (1982) 
friends-of-friends method. Galaxies are linked if sepa- 
rated by a projected distance on the sky, D < Du m , and 
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a velocity difference 5v < 5vu m . In order to minimize 
redshift dependences, we use a sample with a fixed lim- 
iting value M - M* < 2.0. We take 5v Hm = ZbOkms- 1 
and Du m = 520C~ 1 / 2 kpc, where C is the completeness 
of the spectroscopic sample as a function of r magni- 
tude and field position. With these selection parame- 
ters, about half of the galaxies are assigned to groups, 
and half are isolated. For a local comparison, we have 
constructed a group catalog from the NGP sample us- 
ing the same selection criteria (the PM2GC sample is of 
a very narrow strip, making the identification of groups 
problematic) . Although the mass distribution varies less 
with group membership than with local density (Vulcani 
et al. 2012b) we still weight galaxies in the two groups to 
produce similar mass distributions. In contrast to Figure 
7, we see here a clear difference in behavior with environ- 
ment: only isolated galaxies appear to be young, galaxies 
with at least one massive neighbor follow the evolution 
expected for an early formation epoch. 

4. DISCUSSION 

The analysis of the previous section indicates the un- 
expectedly rapid evolution of star formation cannot be 
explained by starbursts, but can be understood as a con- 
sequence of the delayed formation of some fraction of less 
massive galaxies. However, this analysis has not pro- 
vided much quantitative detail on the age distribution 
of galaxies. Given the infinite variety of possible forms 
of star formation histories, the limited information- dis- 
tributions of SSFR's- we have used is insufficient to do 
so. Noeske et al. (2007b) present one possible scenario- 
exponcntially declining Bruzual & Chariot models, like 
those presented in Fig.l, with a delay in formation red- 
shift of the form (1 + Zf orm ) ~ M^ al , but it is hardly 
unique. It is also clearly not correct, since Fig. 6 shows 
that, at every mass, many galaxies appear to be old, and 
only the fraction of younger objects changes with mass, 
and Figure 8 suggests that only isolated galaxies may be 
younger. A more complete description of galaxy forma- 
tion is not possible within the limitations of the data we 
have considered. 

There is, however, more information available. The 
entire 2-dimensional distribution of SSFR-Mass values 
of galaxies over the redshift range 0.0 < z < 1.0, plus 
the integrated star formation history of galaxies back 
the redshifts of order 3, which has been studied by many 
workers (see Cucciati et al. 2012 for references) provides 
additional constraints. In addition, if one makes the com- 
mon assumption that the different structures of bulges 
and disks reflect different timescales for star formation in 
these subsystems, the distribution of bulge-to-disk ratios 
of galaxies provide an additional constraint. In Paper IV 
we will make a start at such a complete analysis: using 
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the distribution of SSFR vs mass and redshift, as well 
as the history of the specific star formation density with 
redshift to constrain the evolution of galaxies with log- 
normal star formation histories. We will, therefore, defer 
all further discussion of this topic to the following paper. 

One item which we will not address in Paper IV is the 
environmental dependance of galaxy ages. Figures 7 and 
8 seem to be telling different stories: galaxy age depends 
on group membership but not on local galaxy density. 
These results are not, however, as inconsistent as they 
appear. Although the two measures of environment are 
correlated, the correlation is very broad: almost half of 
the isolated galaxies fall into the highest two thirds of 
relative density, and 20 percent of the grouped galaxies 
fall into the lowest third of relative density. In general, 
our measure of local density samples a significantly larger 
volume than does our group/isolated criterion. This sug- 
gests that it is the very local environment which affects 
galaxy age: a not implausible conclusion. Together, the 
dependance on galaxy mass and group environment sug- 
gest that galaxy age might be dependent on only a single 
physical parameter: total mass within a volume of size 
about 0.5Mpc. 

The more sophisticated theoretical expectations are 
less than clear. The models of Cen (2011) produce 
more rapid downsizing in denser regions, but Cen does 
not distinguish between archaeological downsizing and 
downsizing in time (to use the terminology of Neistein 
et al. 2006). Neistein et al. _ state that downsizing in 
time, the mass-dependent delay in galaxy formation dis- 
covered by Cowie et al. (1995) cannot be explained with- 
out invoking poorly understood hydrodynamic processes, 
an opinion shared by Fontanot et al. (2009), so the con- 
cordance of theory and observations is not yet estab- 
lished. 
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